Spin Wave Theory of Spin 1/2 XY Model with Ring Exchange on a Triangular Lattice 
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We present the hnear spin wave theory calculation of the superfluid phase of a hard-core boson 
J-K model with nearest neighbour exchange J and four-particle ring-exchange K at half filling on 
the triangular lattice, as well as the phase diagrams of the system at zero and finite temperatures. 
We find that the pure J model {XY model) which has a well known uniform superfluid phase with 
an ordered parameter Mx — (Sf) 7^ at zero temperature is quickly destroyed by the inclusion of 
a negative-K ring-exchange interactions, favouring a state with a (^,0) ordering wavevector. We 
further study the behaviour of the finite-temperature Kosterlitz-Thouless phase transition (Tkt) 
in the uniform superfluid phase, by forcing the universal quantum jump condition on the finite- 
temperature spin wave superfluid density. We find that for i^<0, the phase boundary monotonically 
decreases to T = at K/J = —4/3, where a phase transition is expected and Tkt decreases rapidly 
while for positive K, Tkt reaches a maximum at some K ^ 0. It has been shown on a square lattice 
using quantum Monte Carlo(QMC) simulations that for small K>0 away from the XY point, the 
zero-temperature spin stiffness value of the XY model is decreased^. Our result seems to agree with 
this trend found in QMC simulations. 



I. INTRODUCTION 

The effective studies of continuum field theories have 
resulted in detailed predictions for the low-energy physics 
of quantum spin systems in two dimensions. Spin wave 
theory can provide us with a rather accurate picture of 
the low-lying states of quantum spin systems. There 
are several versions of spin wave theory. The stan- 
dard spin wave theory is based on Holstein-Primakoff 
representation^^ which was first applied to the study of 
Heisenberg model by Anderson^ and further extended to 
second order by Kubo^ and Oguchi^. 

Spin wave theory was thought to be unsatisfactory 
in the case of XY model until Gomez-Santos and 
Joannopoulous^ showed that by a good choice of quan- 
tization axis, one can obtain a good theoretical result 
for XY model. Since then numerous applications of spin 
wave theory have been carried out on XY model with dif- 
ferent lattice configurations and the results obtained so 
far are in a good agree ment with quantum Monte Carlo 
simulations(QMCpII^. 

Another interesting area is the multiple (ring) spin- 
exchange models. This model was first introduced to 
describe the magnetic properties of solid ^He^. It incor- 
porates ring-exchange interactions over plaquettes such 
as spin 1/2 four- spin XY ring exchange of the form: 



Hk = -KY, {StS-StS- + S-Sp-St) , (1) 

{ijkl) 

where the summation runs over plaquettes with the in- 
dices running counter-clockwise, i and j are nearest 
neighbours lying opposite to k and /. 

The ring exchange interaction is important in Wigner 
crystal near the melting densit}'^^. This model, alone or 
in competition with pure XY model with nearest neigh- 
bour exchange has attracted considerable attention over 



the years and has been studied extensively on a square 
lattice using a stochastic series expansion (SSE) quan- 
tum Monte Carlo method^ also a comprehensive theoret- 
ical study (spin wave theory) has been done on a square 
lattice- . It has been suggested in recent works that mod- 
els of this form may harbour exotic ground state prop- 
erties, including de-confined quantum critical points or 
quantum spin liquid phases. 

In this paper, we shall calculate the exact (numeri- 
cal) value of the Kosterlitz-Thouless temperature, super- 
fluid density and other low temperature thermodynami- 
cal properties of spin 1/2 XY model with ring exchange 
interaction (IT]) on a triangular lattice using linear spin 
wave theory without the inclusion of the repulsive inter- 
action between bosons. 

The format of the paper is as follows: In Sec. II, we 
present the model Hamiltonian. In Sec. Ill, we apply lin- 
ear wave theory by choosing our quantization axis along 
the X direction and use it to diagonalize the Hamiltonian 
and obtain its energy. In Sec. IV, we analyze the dis- 
persion and plot it for some values of K/ J . In Sec.V, we 
explore the zero temperature superfluid density using the 
diagonalized Hamiltonian and plot it for some values of 
K/J. In Sec. VI, we calculate the finite temperature su- 
perfluid density and the value of the Kosterlitz-Thouless 
temperature for this model and in Sec. VII, we make 
some concluding remarks. 



II. MODEL 



In this section, we shall present the model Hamiltonian 
and the mean field theory argument of the J-K model on 
a triangular lattice, similar work was done on a square 



lattice^. Our model Hamiltonian is given by 

-K 2_^ [S^ Sj Sj^ Si + 5*^ Sj Sf^ Si ) . 

{ijkl) 



(2) 



The first summation is over nearest neighbour pairs 
on a triangular lattice and the second summation runs 
over the three possible plaquette orientations on a trian- 
gular lattice as shown in Fig.(fT]). It has been shown that 
this Hamiltonian undergoes a Kosterlitz-Thouless phase 
transition for i^ = at Tkt ^0.69 for 2D model and a 
superfluid phase for temperatures less than Ttk^- For 
J > O^K > the Hamiltonian leads to an in-plane quan- 
tum ferromagnet and for K < 0, there is a sign problem 
which prevents (QMC) simulations^ though it is not pos- 
sible to capture this sign problem in the linear spin wave 
theory. QMC simulation on a square lattice has been per- 
formed for i^ > in which there is no sign problem^ ^. 
Lets define the spins as classical vectors by making the 
transformation S~ = pe'^^\ The model Hamiltonian be- 
comes 



H = -2J^p^cos(^i -(l)j) 

-2K ^ p^ cos (0i - (l)j + 0fe - 00 . 

{ijkl) 



(3) 



Now lets consider the case J^K > 0, in this case, the 
minimum energy permits us to have (j)i = (j)j for the J 
term and (j)i = (j)j^ (j)k = 0/ for the K term which leads 
to a ferromagnetic ordering of spins. Consider the case 
J, i^<0, in this case, the minimum energy permits us to 
have (j)i — (j)j = 7T or (j)i = 0, (pj = tt for the J term and 
0i - (/)j + (/)fe - (/>^ = or (t)i = (/)j = (/)fc = 0, (/)^ = TT for 
the K term which leads to ud (up-down) state for the J 
term and uuud state on the plaquettes, basically we have 
two configurations of spins on the lattice. 

This model has also been considered by L. Balents and 
A. ParamekantP^ with the inclusion of the repulsive in- 
teraction U term between bosons and also in the regime 
where J « K. They showed that when J = 0, the 
four-spin exchange leads to a manifold of ground states 
with gapless excitations and critical power-law correla- 
tions and when J 7^ 0, fluctuations select a four- fold 
ferrimagnetically ordered ground state with a small spin 
(superfluid) stiffness which breaks the global U{1) and 
translational symmetry but they did not obtain any exact 
(numerical) value of the Kosterlitz-Thouless temperature 
and the superfluid density. 

From their Hamiltonian they argued that with J = 0, 
the ground state is independent of the sign of K and for 
non-zero J, the sign of K is vital, for J < with i^ < 0, 
there is a ferromagnetic phase while J > leads to a 
^/3 X >/3 Neel order which are also the ground states for 
large \J/K\. They finally concluded that there are no 
phase transitions at any non-zero J other than the well 
known phases. 
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FIG. 1: (Color online) The three plaquette orientations (a), 
(&), and (c) on a triangular lattice with the position coordi- 
nates Sc. = (1,0), s^ = (-h-^)^ ^7 '- 
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III. LINEAR SPIN WAVE THEORY 

The basic assumption of spin wave theory lies on se- 
lecting a classical ground state and determining the fluc- 
tuation around it. In other words, one considers quantum 
fluctuations very close to an ordered ground state con- 
figuration of the system under study. By the usual map- 
ping between spins 5* = 1/2 and the hard-core bosons, we 
can view ([2]) as a hard-core boson model. For J » K 
or i^ = 0, the T = ground state (in-plane ferromag- 
net) has an ordered parameter M^ = (Sf) ^ which 
breaks the U{1) global rotational symmetry or a super- 
fluid phase in the hard-core boson version^^ ^^. One can 
therefore perform a spin wave expansion around this or- 
dered state configuration by introducing the boson oper- 
ators ai and a] which represent the low-energy spin wave 
excitations out of (Sf) and treat other terms in Q as 
perturbations. 

We shall follow the procedure of Gomez-Santos and 
Joannopoulos^ and choose our quantization axis along 
the X direction (instead of z direction). This allows us 
to write the Holstein-Primakoff representatioiP^ in the 
linear spin wave theory as 



Si = 2 -4«i. 



Sf 



1 f t 

— aj — ai 
2i\ ' 



(4) 




The method for diagonalizing the Hamiltonian |2| 
given by R. Schaffer et al^ . It involves writing (p 
terms of S^ and 5'J using S^ = S^ ^ iSj^ after which 
we substitute Q into the resulting equation and then 
Fourier transform taking into account the summation 



over all the three plaquettes for the i^-term and finally 
diagonalize using Bogoliubov transformation. After all 
these mathematical gymnastics, the diagonalized Hamil- 
tonian is 



H = Hmf + ^ (cJk - Ak) 

k 



(5) 



Here, the mean-field energy and the coefficients are 
totally different from those obtained on a square lattice^, 
they are given by 



H^ 



MF 



-3(-JiV- 



1 



KN 



Ak = JQk + KRk, 
Bk = JSk + KTy,, 



where 



Qv 



3 1 



7k^ 
2 

Jk = X7k, 



(-1) 



i?k = 3 <( ^ - -7k + - (7k + 7k) 



Tv 



= 3 I -7k - g (7k + 7k) 



(6) 

(7) 
(8) 
(9) 



(10) 

(11) 
(12) 

(13) 



and the lattice structure constants are given by 
1 



7k 



7k 



cos kr 



2 cos — cos 

2 2 



1 / /-^ 3kx VSky \ 
- cos v6ky -\- 2 cos cos 



(14) 



We can see that the Hamiltonian ([5| reduces to pure XY 
model in the limit iC = as expected. This diagonalized 
Hamiltonian ([5| will be used to analyse some properties 
such as dispersion, ground state, internal energy, super- 
ffuid density etc., all as a function of Kj J . 



IV. DISPERSION 

We shall now study the dispersion of J-K model as a 
function of Kj J . Similar to J model, the dispersion is 
given by 



e(k)=2c^k = 2jA2-52 



k' 



(15) 



but in this case the coefficients A\^ and B\^ are given by 
(|8| and (|9|. The graph of ([l5| is shown in Fig.(|3| for 



specific values of Kj J . For Kj J = 
takes the compact form 



-4/3, the dispersion 



e(k) = vr^^, 



(16) 

where 7k is given by ( p3| ). For K/J = — 2, the dispersion 
e (k) = a/O = for all values of k and no reasonable plot 
can be obtained below this value of K/J because the 
dispersion develops an imaginary part along the kx and 
ky directions. The dispersion is also zero at k = (0, 0) 
for ah values of K/J. At k = (^,0),(^,^) which 
corresponds to the corners of Brillouin zone Fig.Q, it is 
given by 



e{k) = -y/24- 



m{K/J) + [ZK/Jf 



(17) 



which clearly vanishes at two points K/J = —4/3 and 
K/J = -2. The plot of ^ is shown in Fig.Q. There 
are two regions I and II, region I corresponds to K/J > 
—4/3 and region II corresponds to K/J < —2. Since the 
dispersion ( 15 ) develops an imaginary part along the kx 
and ky directions for K/J < —2, we see that region II is 
not physical. 




FIG. 2: (Color online) The dispersion e (k) as 
K/J at the corners of the Brillouin zone 



function of 



The dispersion along the kx direction is also shown in 
Fig. ([3]) for several values of K/J. In this case the dis- 
persion shows a linear function of kx and a zero mode 
at the center of the Brillouin zone k=Q=(0,0) except 
for K/J = —4/3 where it develops three zero modes 
(minima) at k= ±Q =(±^,0) and k= Q=(0,0). If 
one chooses ferromagnetic ordering along the ky axis 
then the corresponding ordering wave vector is k= ±Q 
= (±^,0) which is the soft modes of the dispersion for 
K/J = —4/3. The linear spin wave instability of the ex- 
citation spectrum at the corner of the Brillioun zone k= 
±Q =(±^,0) occurs for K/J = -4/3. For the pure 
XY model K/j = 0, there is a gapless excitation at Q 
= (0,0) (Goldstone mode of the superfiuid phase), but 
there is no minima at the ordering wave vector. Near the 
zero modes for K/J = —4/3 we have 



e(k) 



V^ 



[k^ + kt 



.1/2 



VSJk. 



(18) 






FIG. 3: (Color online) The dispersion e (k) along kx direc- 
tion for several values of K/J (first), the dispersion e (k) as a 
function of k = {kx, ky), for K/J = 4/3 (second) and K/J = 
-4/3 (third). 
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FIG. 4: (Color online) The Brillouin zone of the triangular 
lattice. 



using the labeling in Fig.Q, we have Oi = Oi^ Oj = Ok-, 
therefore the dependence of the ring exchange terms 
cancel. The twisted J-K Hamiltonian becomes 

H{e) = -2 J ^ { (^f ^J + S^S]) COS 6 

{ij) 

-^{SfS] -S^S^) sin 0} (20) 

{ijkl) 



In the linear spin wave theory ( |2Q| ) reduces to 
H{0) = -2J^ {SfS^ + S^S^) cosO 

{ij) 

-kY, {sfs-sts^ + s-sp^s^) . 

{ijkl) 



(21) 



Notice that the coefficients of sin 6 vanishes in the case. 
Comparing ([2| and (21 ), we see that the effect of the twist 
is to rescale the nearest-neighbour exchange interaction 
inpl) as J ^ J cos 6. Therefore, the diagonalized form of 
(|21|^is exactly of the form ([5| with J replaced by J cos 6. 
That is 



H{0) = Hmf{0) + Yl (^k(^) - A^{0)) 

k 

+ Y ^k(^) (^<^k<^k + Q^Lk<^-k j , 



(22) 



V. SUPERFLUID DENSITY AT T = 

In this section we shall explore the superfiuid density 
of the full J-K model on a triangular lattice. In this case 
we must consider the influence of the ring exchange term 
on the applied phase twist 0. The relevant transforma- 
tions are similar to that of a square lattice^^, with this 
transformations the ring exchange term becomes 

Sp-Spf = 5+^-5+5^6^^^^+^^-^^-^^), (19) 



where the twisted mean-field energy and the coefficients 
are given by 



Hmf{0) = -3( ^JNcosO^^KN 



Ak(6>) = JQkCosi9 + Ki?k, 
B^,{e) = JS^,cosO^KT^,, 



(23) 
(24) 
(25) 



the coefficients Qic^k^^'k, and Tk remain the same as 
([llj-([l3|. 

The free energy is given by 



-^ek(^) 



where the ground state energy Eo{0) is given 

Eo{e) = HMFie) + J2 (^k(^) - A^{e)) . 



(26) 



(27) 



At zero temperature j3 = oo, the free energy is simply the 
ground state energy. We shall calculate the superfluid 
density from the first principle. Taylor expanding the 
ground state energy we have 



^o(^) = Eo{e = 0) + ^p. 



■O{0% 



Ps{T = 0) 



1 d^Eoje) 
N 56*2 ' 



(28) 



(29) 



Prom (27) and (29) we obtain 
J 



P«(T = 0) = 



3J 



2N 



{J{Ql 



si 



E[Q.-- (30) 

K{QuRk-S^Tu)^ 



Note that we have divided by 2 to account for the 
dimensionahty of the lattice. 

The plot of ps{T = 0) is shown in Fig.([5]) for a range 
of K/J^ we have set J = 1/2. The superfluid density 
curve has its maximum at K/J = 0, with a value of 
ps{T = 0) = 0.4059 and decreases monotonically with 
increasing \K\ as one moves away from this maximum, 
in other words the ring exchange term decreases the value 
of ps{T = 0) from the pure XY result. On the positive-iC 
side, ps decreases relatively gradually, only becomes zero 
for extremely large values. This is consistent with the 
result obtained in the dispersion, which indicates that no 
soft modes develop for moderate values of positive K. 
Similar result was observed on a square lattice^, which 
is consistent with quantum Monte Carlo simulation. On 
the negative- i^ side, the value of ps decreases rapidly as 
it approches K/J = —4/3. 



VI. SUPERFLUID DENSITY AND 

KOSTERLITZ-THOULESS TRANSITION AT 

FINITE TEMPERATURE 

In this section we shall discuss the uniform superfluid 
phase in the J — K model at finite temperatures. In order 
to calculate the finite temperature superfluid density, we 
replace the twisted ground state energy in (29) with the 
twisted free energy, that is 



Ps{T^O) 



1 d'^F{0 ) , 
N~dW 



\0=O- 



(31) 
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FIG. 5: (Color online) The superfluid density as a function of 
K/J. The dashed line is the mean- field result, p^^ = 0.375. 
The linear SW result for K/J = is pf ^(T = 0) = 0.4059. 



Using (26) we obtain 



PsiT^O) 



3J 

T 



j_ 

2iV 



e[q^ 



1 



(32) 



(^kQk - 5k^k) 1 



oek/T _ I 



This expression for psiT) is plotted as a function of T / J 
for several values of K/J (see Fig.(|6])). 
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FIG. 6: (Color online) The superfluid density as a function 
of T/J. Lines are labelled by the parameter value K/J. The 
dashed line is the universal jump condition. 

The graph shown in Fig.(|6| is similar to the one ob- 
tained on a square lattice^ which shows slowly decay- 
ing superfluid density. The dash line is the so-called 
universal jump condition 






(33) 



which accounts for the discontinuity in ps{T). The esti- 
mate of Tkt can be found by solving ps(T)/T = 2/7r, 
using PsiT) from our spin wave theory. Using this 
procedure we found the KT transition temperature at 
Tkt = 0.9295J for K = and J = 1/2, the parameter 
values for the pure XY model. This is shown in Fig. (pi). 




FIG. 7: (Color online) The KT transition phase boundary as 
calculated from linear SW theory. A maximum of Tkt/J = 
1.201 occurs at approximately K/J = 6.019. 



where the dashed line crosses the curve for K/J = 0. 
The plot of Tkt for non-zero values of K/J is shown in 
Fig.([7|). It is interesting to see that the maximum does 
not occur at K/J = but at K/J ^ 6.019 before drop- 
ping slowly. Similar to the case of square lattice^, the 
KT transition drops to zero for extremely large value of 
K/J which cannot be captured by a simple spin wave the- 
ory. On the negative- i^ side, the phase boundary drops 
rapidly as K/J approaches —4/3. 



VII. CONCLUSION 



this model was calculated by applying the traditional 
Holstein-Primakoff representation and summing over the 
three plaquettes orientations on a triangular lattice. This 
calculation showed a spin wave instability and a develop- 
ment of three minima at k= ±Q = (± ^ , O) , and k = Q= 
(0, 0) for K/J = —4/3. One should expect a phase tran- 
sition from a superfluid phase to another phase at this 
wave vector k= Q =(^,0) for K/J = —4/3. A more 
careful analysis of quantum Monte Carlo data should pro- 
vide further insight into this issue. 

The mean field superfluid density and ground state 
(zero temperature) spin wave superfluid density obtained 
in this model for K/J = is bigger than that of a square 
lattice^ despite the fact that both are 2D systems. This 
might be due to a larger number of nearest neighbours 
and plaquettes on a triangular lattice. We calculated the 
finite temperature uniform superfluid density and use it 
to estimate the Kosterlitz-Thouless transition tempera- 
ture by forcing it to obey the universal quantum jump 
condition. We find that for i^<0, the phase boundary 
monotonically decreases to T = at K/J = —4/3, where 
a phase transition is expected. It has been shown with 
the Model Hamiltonian ( 2]) on a square lattice using QMC 
simulations that for small K>0 away from the XY point, 
the zero-temperature spin stiffness value of the XY model 
is decreasecP. Our results above seem to agree with this 
trend found in QMC simulations. 



We have presented a comprehensive study of linear spin 
wave theory of hard-core bosons (zero field XY model) at 
half filling on the triangular lattice. We studied through 
linear spin wave theory the destruction of uniform su- 
perfluid phase in the bosonic ring exchange model on 
a triangular lattice at half filling. The dispersion of 
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